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Dynamical and transport properties of a simple single-band spin-fermion lattice model for (III,Mn)V diluted 
magnetic semiconductors (DMS) is here discussed using Monte Carlo simulations. This effort is a continuation 
of previous work (G. Alvarez et al, Phys. Rev Lett. 89, 277202 (2002)) where the static properties of the 
model were studied. The present results support the view that the relevant regime of J/t (standard notation) 
is that of intermediate coupling, where carriers are only partially trapped near Mn spins, and locally ordered 
regions (clusters) are present above the Curie temperature Xc. This conclusion is based on the calculation of 
the resistivity vs. temperature, that shows a soft metal to insulator transition near Xc, as well on the analysis of 
the density-of- states and optical conductivity. In addition, in the clustered regime a large magnetoresistance is 
observed in simulations. Formal analogies between DMS and manganites are also discussed. 
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I. INTRODUCTION 

Diluted magnetic semiconductors (DMS) are attracting 
much attention lately due to their potential for device applica- 
tions in the growing field of spintronics. In particular, a large 
number of DMS studies have focused on III-V compounds 
where Mn doping in InAs and GaAs has been achieved us- 
ing molecular beam epitaxy (MBE) techniques. The main 
result of recent experimental efforts is the discoverjiUAi 
of ferromagnetism at a Curie temperature Xc^llO K in 
Gai-^Mn^As, with Mn concentrations x up to 10%. It is 
widely believed that this ferromagnetism is "carried induced", 
with holes donated by Mn ions mediating a ferromagnetic 
interaction between the randomly localized Mn 2+ spins. In 
practice, anti-site defects reduce the number of holes n from 
its ideal value n=x, leading to a ratio p=(n/x) substantially 
smaller than 1 . 

Until recently, theoretical descriptions of DMS materials 
could be roughly classified in two categories. On one hand, 
the multi-band nature of the problem is emphasized as a cru- 
cial aspect to quantitatively understand these materials 
In this context the lattice does not play a key role and a con- 
tinuum formulation is sufficient. The influence of disorder 
is considered on average. On the other hand, formulations 
based on the possible strong localization of carriers at the Mn- 
spin sites have also been proposed^iS In this context a single 
impurity-band description is considered sufficient for these 
materials. Still within the single-band framework, but with 
carrier hopping not restricted to the Mn locations, recent ap- 
proaches to the problem have used dynamical mean-fieldi^ 
or reduced-basis 13 approximations. An effective Hamiltonian 
for Gai ^Mn^As was derived in the dilute limit and studied in 
Ref . [3 All these calculations are important in our collective 
effort to understand DMS materials. 

However, it is desirable to obtain a more general view of 
the problem of ferromagnetism induced by a diluted set of 
spins and holes. To reach this goal it would be better to use 
techniques that do not rely on mean-field approximations and, 



in addition, select a model that has both the continuum and 
impurity-band formulations as limiting cases. Such an ap- 
proach would provide information on potential procedures to 
further enhance Xc and clarify the role of the many parame- 
ters in the problem. In addition, these general considerations 
will be useful beyond the specific details of Gai-^Mn^As, al- 
lowing us to reach conclusions for other DMS. An effort in 
this direction was recently initiated by Mayr and two of the 
authors. 15 A detailed Monte Carlo study of a simple model 
for DMS already allowed us to argue that Xc could be further 
enhanced -perhaps even to room temperature- by further in- 
creasing x and p from current values in DMS samples. These 
predictions seem in agreement with recent experimental de- 
velopments since very recently Gai-^Mn^As samples with 
Xc as high as 150 K were prepared^ a result believed to be 
caused by an enhanced free-hole density. Also samples with 
Xc ~ 127 K and x > 0.08 were reported in Ref. 17, and, very 
recently, a Xc of 140K was achieved on high quality GaMnAs 
films grown with Arsenic dimers; 18 There seems to be plenty 
of room to further increase the critical temperatures according 
to these theoretical calculations. 

The model used in Ref.fl5lwas a single-band lattice Hamil- 
tonian, also studied by other groups, which does not have 
the multi-band characteristics and spin-orbit couplings of the 
real problem. However, it contains spin and holes in inter- 
action and it is expected to capture the main qualitative as- 
pects of carrier-induced ferromagnetism in DMS materials. 
The choice of a single-band model allows us to focus on the 
essential aspects of the problem, leaving aside the numeri- 
cal complexity of the multi-band Hamiltonian, which unfor- 
tunately cannot be studied with reliable numerical techniques 
at present without introducing further approximations. 19 Note 
that a quite similar philosophy has been followed in the re- 
lated area of manganites for many years. In fact, it is well- 
known that the orbital order present in those compounds can 
only be studied with a two-band model. However, single-band 
approaches are sufficient to investigate the competition be- 
tween ferromagnetic and anti-ferromagnetic states, an effect 
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FIG. 1: Phase diagram of the single-band model Eq. Q, as discussed 
in Ref.[l3 The figure shows schematically the three coupling regions 
discussed in the text: (I) weak, (II) intermediate, and (III) strong cou- 
pling, as well as the Tc and T* dependence on the coupling J/t. The 
region between T* and Tc contains FM "clusters", with magnetic 
moments that are not aligned. 



recently argued to be at the heart of the colossal magnetore- 
sistance phenomenon As a consequence, a qualitatively 
reliable study of the single-band DMS model - with itinerant 
holes and localized classical spins - is expected to be impor- 
tant for progress in the DMS context as well. 

In our previous publication^ the phase diagram of the 
single-band model for DMS was already sketched, 22 This 
phase diagram is reproduced in Fig [I] The only parameter 
in the Hamiltonian is the ratio J/t, where J is the local anti- 
ferromagnetic coupling between the spins of carriers and lo- 
calized Mn 2+ ions, and t is the carrier-hopping amplitude (for 
details, see the full form of the Hamiltonian in the next sec- 
tion). Due to its large value 5/2, the Mn spin is assumed 
to be classical for numerical simplicity. At small J/t, in- 
dividual carriers are only weakly bounded to the Mn spins. 
Even for small realistic values of x, the carrier wave func- 
tions strongly overlap and a large-bandwidth itinerant band 
dominates the physics. In this regime, Tq grows with J/t. 
This portion of the phase diagram is referred to as 'weak cou- 
pling' in Fig[T] and in our opinion it qualitatively c orre sponds 
to the continuum-limit approach pursued in Refs. l5tf5l7t In the 
other extreme of large J/t, the previous effort 15 found a Tc 
that is much suppressed, converging to zero as J / 1 increases 
(at least at small values of x). The reason for this behavior 
is the strong localization of carriers at Mn spin locations, to 
take advantage of the strong J coupling. The localization sup- 
presses the mobility and the carrier-induced mechanism is no 
longer operative. However, when several Mn spins are close 
to one another small regions can be magnetized efficiently. As 
a consequence, a picture emerges where small islands of ferro- 
magnetism are produced at a fairly large temperature scale T* 
(>Xc) that grows with J, but a global ferromagnetic state is 
only achieved upon further reducing the temperature such that 
the overlaps of wave functions induce a percolation-like pro- 
cess that aligns the individual preformed clusters The proper 
procedure to distinguish between these two regimes is through 
an analysis of short- and long -dista nce spin correlations. In 
our opinion, the effort of Refs. l9ll0l belongs to this class, and 



it corresponds to the 'strong coupling' region of Fig [I] Be- 
tween the weak and strong coupling domains, an intermedi- 
ate region provides a natural interpolation between those two 
extreme cases. The Curie temperature is optimal (i.e., the 
largest) in this intermediate zone, at fixed values of p and x. In 
this regime, there is sufficient interaction between the clusters 
to become globally ferromagnetic, and at the same time the J 
is sufficiently strong to induce a robust Tc . If chemical con- 
trol over J/t were possible, this intermediate region would be 
the most promising to increase the Curie temperature (again, 
at fixed x andp). 

However, it may even occur that the DMS materials are 
already in the intermediate optimal range of J/t couplings. 
An indirect way to test this hypothesis relies on calculations 
of dynamical and transport properties, which are presented 
in this paper, and its comparison with experiments. In par- 
ticular, the experimentally measured d.c. resistivity pdc of 
DMS compounds has a nontrivial shape with a (poor) metal- 
lic (dpdc/dT>0) behavior below Tc, which turns to insu- 
lating (dpd c /dT<0) at higher temperatures. This nontrivial 
profile resembles results reported in the area of manganites, 
which also have a metal-insulator transition as a function of 
temperature, 21 although in those materials the changes in re- 
sistivity with temperature are far more dramatic than in DMS. 
The formal similarities DMS-manganites have already been 
remarked in previous literature^ and suggest a common ori- 
gin of the pdc vs. temperature curves. In Mn-oxides it is 
believed that above the Curie temperature, preformed ferro- 
magnetic clusters with random orientations contribute to the 
insulating behavior of those materials For DMS, a simi- 
lar rationalization can be envisioned if the state of relevance 
above Tc has preformed magnetic clusters. From the previous 
discussion 15 it is known that the intermediate J/t region can 
have clusters without collapsing Tc to a very small value. At 
the Curie temperature, the alignment of these preformed mo- 
ments leads to a metallic state. An explicit calculation of the 
resistivity - using techniques borrowed from the mesoscopic 
physics context - is provided in this paper, and the results sup- 
port the conjecture that clustered states can explain the trans- 
port properties of DMS. These unveiled analogies with man- 
ganites are not just accidental. Clearly, approaches to DMS 
that rely on mobile carriers and localized spins in interaction 
have close similarities with the standard single-band double- 
exchange model where only one e 9 -orbital is considered. The 
key difference is the presence of strong dilution in DMS as 
compared with manganites. 

It is interesting to remark that in agreement with the 'clus- 
tered' state described here, recently Timm and von Oppen 
have shown that correlated defects in DMS are needed to de- 
scribe experimental data. 23 Their simulations with Coulombic 
effects incorporated lead to cluster formation, with sizes well 
beyond those obtained from a random distribution of Mn sites 
as considered here. In this respect, the results of Ref.|23| pro- 
vide an even more dramatic clustered state than ours. If the 
Mn spins were not distributed randomly in our simulations but 
in a correlated manner, the state above Tc would be even more 
insulating than reported below and the physics would resem- 
ble much closer that observed in manganites. Note also that 
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Kaminski and Das Sarma have also independently arrived to 
a polaronic state 24 that qualitatively resembles the clustered 
state discussed here. It is also interesting that in recent ion- 
implanted (Ga,Mn)P:C experiments that reported a high Curie 
temperature, the presence of ferromagnetic clusters were ob- 
served and they were attributed to disorder and the proximity 
to a metal-insulator transition. 25 Even directly in (Ga,Mn)As, 
the inverse magnetic susceptibility deviates from the Curie 
law at temperatures above Tq (Fig. 2(b) of Ref.0), which may 
be an indication of an anomalous behavior. In manganites, 
these same anomalies in the susceptibility are indeed taken as 
evidence of the formation of clusters^ 

The paper is organized as follows: In Section[n]the details 
of the method of simulation, definitions, and conventions are 
presented. Section |In| describes the results for the density - 
of-states (DOS) of this model, the appearance of an "impu- 
rity band" at large J/t, and the relation of the DOS with 
the optimal Tq. The temperature and carrier dependence of 
the optical conductivity and Drude weight are presented in 
Section |lVl with results compared with experiments. In Sec- 
tion 13 the resistance of a small cluster is calculated and also 
compared with experimental data. The intermediate coupling 
regime is the only one found to qualitatively reproduce the 
available resistivity vs. temperature curves. From this per- 
spective, both the weak and strong coupling limits are not re- 
alistic to address DMS materials. In Section EU the influence 
of a magnetic field is studied by calculating the magnetore- 
sistance. Large effects are observed. Section IVlTl describes 
the characteristics of the clustered- state regime above Tq at 
strong coupling. Finally, in Section lVIlTl the 'transport' phase 
diagram is presented, and it is discussed in what region it ap- 
pears to be the most relevant to qualitatively reproduce results 
for Gai-^Mn^As MBE-grown films. 



II. METHOD AND DEFINITIONS 
A. Hamiltonian Formalism 

There are two degrees of freedom in the DMS problem de- 
scribed here: (i) the local magnetic moments corresponding to 
the 5 electrons in the d shell of each Mn impurity, with a total 
spin 5/2, and (ii) the itinerant carriers produced by the Mn im- 
purities. Since these Mn impurities substitute Ga atoms in the 
zinc-blende structure of GaAs, there is, in principle, one hole 
per Mn. However, it has been found that the system is heavily 
compensated and, as a consequence, the actual concentration 
of carriers is lower. In the present study, both the density of 
Mn atoms, x, and the density of carriers, n, are treated as in- 
dependent input parameters. Moreover, the ratio p = n/x is 
defined, which is a measure of the compensation of the sys- 
tem, e.g. for p = the system is totally compensated and for 
p = 1 there is no compensation. 

The Hamiltonian of the system in the one-band approxima- 
tion can be written as: 

H = -t £ c^ + J^Srf;, (1) 

<ij>,a I 



where c] a creates a carrier at site i with spin a. The carrier- 
spin operator interacting antiferromagnetically with the lo- 
calized Mn-spin Si is 07 = c\ a a a ^ci/3. Through nearest- 
neighbor hopping, the carriers can hop to any site of the 
square or cubic lattice. The interaction term is restricted to 
a randomly selected but fixed set of sites, denoted by /. Note 
that in addition to the use of a single band, there are other 
approximation in the model described here: (i) The Hubbard 
U/t is not included. This is justified based on the low-carrier 
concentration of the problem, since in this case the proba- 
bility of double occupancy is small. In addition, Mn-oxide 
investigations 2 ^ 2 ! have shown that an intermediate or large J 
coupling acts similarly as U/t, also suppressing double occu- 
pancy, (ii) A nearest-neighbors anti-ferromagnetic coupling 
between the Mn spins is not included. Again, this is justified 
as long as x is small, (iii) Finally, also potential disorder is ne- 
glected (together with the spin, the Mn sites should in princi- 
ple act as charge trapping centers). This approximation is not 
problematic when the Fermi energy is larger than the width of 
the disorder potential (i.e. large p and x). However, in the (p, 
x) range analyzed in this paper it becomes more questionable 
to neglect this effect. Its influence will be studied in a future 
publication. 

In the present study the carriers are considered to be elec- 
trons and the kinetic term describes a conduction band. How- 
ever, in Gai-^Mn^As the carriers are holes and the valence 
band has to be considered instead. Although the latter is more 
complicated, the simplified treatment followed here should 
yield similar results for both cases £A2iI£ 

The local spins are assumed to be classical which allows 
the parametrization of each local spin in terms of spherical 
coordinates: (Oi,<f>i). The exact-diagonalization method for 
the fermionic sector is described in this section, following 
Ref.|2(]. The partition function can be written as: 

z = n d9i sin ^ I n & ») • (2) 

where Z g ({0i,(j)i}) = Tr(e -/3K ) is the partition function of 
the fermionic sector, K = H — fiN, N is the number operator, 
and /i is the chemical potential. 

In the following, an hypercubic lattice of dimension D, 
length L, and number of sites N = L D will be considered. 
Since K is a Hermitian operator, it can be represented in terms 
of a hermitian matrix which can be diagonalized by an unitary 
matrix U such that 
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The basis in which the matrix K is diagonal is given by the 
eigenvectors i4|0), • • •> u \n\®)> where the fermionic opera- 
tors used in this basis are obtained from the original operators 
through u m =J2j a U^jpCja, with m running from 1 to 2 N. 
Defining ^J n ^ m =n m and denoting by n m the eigenvalues 
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of n m , the trace can be written 



Tr g (e (3K ) = ^2 (ni---n 2N 



-(3K\ 



ni,...,ri2N 



= Yl {ni..-n2N\e-^=^\n l ...n 2 N), (4) 



ni,...,n 2 N 



since in the {i4n|0)l basis, the operator K is Yl\ e A^A> an d 
the number operator can be replaced by its eigenvalues. Thus, 
Eq. can be rewritten as 

Z = Y[( MiShiOi d&)T[(l + e-^), (5) 

• Jo Jo x=1 

which is the formula used in the simulations. 



B. Monte Carlo Method 

The integral over the angular variables in Eq. © can be per- 
formed using a classical Monte Carlo simulation^ The eigen- 
values must be obtained for each classical spin configuration 
using library subroutines. Finding the eigenvalues is the most 
time consuming part of the numerical simulation. The inte- 
grand is clearly positive. Thus, "sign problems" in which the 
integrand of the multiple integral under consideration can be 
non-positive, are fortunately not present in our study. 

Although the formalism is in the grand-canonical ensem- 
ble, the chemical potential was adjusted to give the desired 
carrier density, n. To do so, the equation — n = was 
solved for fi at every Monte Carlo step by using the Newton- 
Raphson method. 27 This technique proved very efficient in ad- 
justing with precision the desired electronic density and, as a 
consequence, our analysis can be considered as in the canoni- 
cal ensemble as well, with a fixed number of carriers. 

Usually, 2,000-5,000 Monte Carlo iterations were used to 
let the system thermalize, and then 5,000 -10,000 additional 
steps were carried out to calculate observables, measuring ev- 
ery five of these steps to reduce autocorrelations. 

Since for fixed parameters there are many possibilities for 
the random location of Mn impurities, results are averaged 
over several of these disorder configurations. Approximately 
10-20 disorder configurations were generally used for small 
lattices (4 3 and 10x10) and 4-8 for larger lattices (6 3 and 
12x12). The inevitable uncertainties arising from the use of 
a small number of disorder configuration does not seem to 
affect in any dramatic way our conclusions below. This is de- 
duced from the analysis of results for individual disorder sam- 
ples. The qualitative trends emphasized in the present paper 
are present in all of these configurations. 



C. Observables 

Quantities that depend on the Mn degrees of freedom (Oi 
and <j>i in the previous formalism) are calculated simply by 
averaging over the Monte Carlo configurations. Note that 



any observable that does not have the continuous symmetry 
of the Hamiltonian will vanish over very long runs. Thus, it 
is standard in this context to calculate the absolute value of 



the magnetization, \M\ = y Si • Sj, as opposed to the 
magnetization vector. Another useful quantity is the spin- spin 
correlation, defined by: 
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(6) 



where N(x) is the number of non-zero terms in the sum. The 
correlation at a distance d is averaged over all lattice points 
that are separated by that distance, but since the system is di- 
luted, the quantity must be normalized to the number of pairs 
of spins separated by d, to compare the results for different 
distances. 

The observables that directly depend on the electronic de- 
grees of freedom can be expressed in terms of the eigenvalues 
and eigenvectors of the Hamiltonian matrix K. The DOS, 
N{u), is simply given by ^2 x S(co — e\). However, the ma- 
jority and minority DOS, N^(lu) and Ni(w), were also cal- 
culated in this study. N^(u) indicates the component that 
aligns with the local spin, i.e., N^(u) is the Fourier trans- 
form of < 4r(£)c;j(0) >, where = cos(0;/2)c^ + 
sin(0 i /2)e- icf)i c ii . Then: 
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N,(w) = ^<5( W -e A )[^/7^ A /7 A , iT cos 2 (^/2) + 

A i 

Ul tX Ux,n sin 2 (fV2) + 

( U h,\ U Kil exp(-^i) + U} ux U X tf exp(z<&)) x 
cos(0i/2)sin(0i/2)], (7) 

where for sites i without an impurity Oi = fa = is assumed. 
A similar expression is valid for Ni(u). The optical conduc- 
tivity was calculated as: 
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where the current operator is: 



<J x (t)-j x (0)>, (8) 



3x = ^Zjfi)^,a c 3^ ~ H - C ')i 



(9) 



with x the unit vector along the x-direction. For uo 7^ 0, cf(uj) 
can be written as: 



(J(CJ) 
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(10) 



(1 _|_ 6 £Oa-aO)(i -f 6 -/3(p a /-a0) 
5(u + px-p\>). (11) 

Both N(u) and a(co) were broadened using a Lorentzian 
function as a substitute to the (5-functions that appear in 



5 



Impurity 
band 



D.O.S. 



2J 




Main band 



Energy 



FIG. 2: Schematic representation of the DOS, for a ferromagnetic 
configuration and strong enough J/t coupling. The "impurity band" 
has weight x, and the chemical potential \i lies within it. 



Eqs. Q and iTTl . The width of the Lorentzian used was 
e = 0.05 in units of the hopping, t. 

The optical conductivity in d dimensions obeys the sum 
rule: 



D _ Tre 2 < -f > 
~2 2Nd 



/ *M 

J0+ 



duo, (12) 



where D is the Drude weight and T is the kinetic energy: 

-T = *E (clcjv + H.c). (13) 



<ij>,a 



Although the Drude weight gives a measure of the d.c. con- 
ductivity properties of the cluster, the "mesoscopic" conduc- 
tance was also calculated in order to gather additional infor- 
mation. The details are explained in Section fvl 



III. DENSITY OF STATES 

In this section, the density-of- states calculated at several 
couplings is shown and discussed. A few basic facts about 
the DOS of the model, Eq. |T}, are presented first. If carrier 
localization effects are not taken into account and a ferromag- 
netic state is considered at strong J/t coupling, then two spin- 
split bands ("impurity bands") will appear for this model, each 
with weight proportional to x, at each side of the unperturbed 
band which would have weight 2(1 — x) (Fig.|2|). For partially 
compensated samples, p < x, the chemical potential will be 
located somewhere in the first spin-split band, and, as a conse- 
quence, only this band would be relevant in the model Eq. |T}. 

In practice, Monte Carlo simulations show that, in the 
regime of interest where Tq is optimal, the impurity band is 
not completely separated from the unperturbed band. For this 
reason, both the unperturbed band and the impurity band have 
to be considered in quantitative calculations. To illustrate this, 



Figure [3] shows N(oj) vs. uj for different J/t's in two dimen- 
sions. With growing J/t, the "impurity band" begins to form 
for J/t > 3.0, and at J/t = 6.0 it is already fairly sepa- 
rated from the main band. However, note that for this extreme 
regime of J / t, Tn is far lower than for J/t = 3.0, as discussed 
before in Ref.fl51 In fact, the optimal J/t for the (x,p) pa- 
rameters used here was found to be in the range 2.0-4.0. This 
implies that the model shows the highest Tq when the "impu- 
rity band" is about to form, but it is not yet separated from the 
main band. This introduces an important difference between 
our approach to DMS materials, and the theoretical calcula- 
tions presented in Ref.^J. Figure|4ja)-(c) indicates that a sim- 
ilar behavior is found in three dimensions. In this case, once 
again, the optimal Tq occurs for intermediate J/t, as seen in 
Fig-Etd) where the phase diagram obtained from a 6 3 lattice 
with x=0.25 andp=0.3 is shown. 

The physical reason for this behavior is that at very large 
J/t , the states of the impurity band are highly localized where 
the Mn spins are. Local ferromagnetism can be easily formed, 
but global ferromagnetism is suppressed by the concomitant 
weak coupling between magnetized clusters. 

Both, Figures [3 and |4| represent results obtained at partic- 
ular values of parameters x=0.25, p=0.3, and T/t = 0.01. 
However, results for several other sets {(x,p, T)} were gath- 
ered as well (not shown). The x dependence of the DOS is 
simple: increasing x produces a proportional increase in the 
number of states of the impurity band, and a corresponding 
decrease of the main band weight (see Fig|2|). In terms of p, in- 
creasing the effective carrier concentration by a small amount 
was found to simply shift the chemical potential to the right. 
Concerning the temperature dependence of the DOS, at low 
temperature the states contributing to the impurity band are 
polarized, i.e. the system is ferromagnetic, as shown by the 
different weights of the majority and minority bands (Figs. [3] 
and|4j. As the temperature increases, spin disorder grows due 
to thermal fluctuations, and when T ~ T*, i.e. in the param- 
agnetic regime, the bands become symmetric. 

It appears that the only change that a small increase in the 
carrier concentration produces is an increase in the chemi- 
cal potential. However, experimentally the appearance of a 
pseudo-gap at the Fermi energy has been reported. 29 Although 
the finite size effects in the present theoretical calculation do 
not allow for a detailed study of N(uj) at or very near the 
Fermi energy with high enough precision, in certain cases a 
pseudo-gap was indeed observed in the present study at the 
chemical potential, particularly when the system is in the clus- 
tered regime. While this result certainly needs confirmation, it 
is tempting to draw analogies with the more detailed calcula- 
tions carried out in the context of manganites, where the pres- 
ence of a pseudo-gap both in theoretical models and in exper- 
iments is well established^^iiSi Pseudo-gaps are also present 
in under-doped high temperature superconductors. Given the 
analogies between DMS materials and transition-metal oxides 
unveiled in previous investigations p£ it would not be surpris- 
ing that DMS presents a pseudo-gap in the clustered regime 
as well. More work is needed to confirm these speculations. 

Results presented in this section for the DOS qualitatively 
agree with those found in Ref . [l2| using the dynamical mean- 
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field technique (DMFT) in the coupling regime studied there. 
However, the DMFT approach is a mean-field approximation 
local in space and, as a consequence, the state emphasized 
in this paper -with randomly distributed clusters- cannot be 
studied accurately with such a technique. 

Experimental evidence of the formation of the impurity 
band for Gai T M n. T As has been provided by Okabayashi et 
al. in Refs. l3lll32n see also Ref.l33l). It is interesting to remark 
that half the total width of the impurity band obtained in our 
simulations is about 2t (see, e.g., Figs. |3]and|4|), and using 
£=0.3eV 15 this is equal to 0.6eV, in good agreement with the 
width 0.5eV estimated from the photoemission experimental 
measurements mentioned before ^ 

The local DOS is studied next. The purpose of this analysis 
is to further understand the inhomogeneous state that forms 
as a consequence of Mn-spin dilution and concomitant carrier 
localization. The local DOS is shown in Fig. |5]on an 8x8 
lattice at x = 0.25, p = 0.4, J/t = 2.5, and T/t = 0.05. In 
this particular case only, an arrangement of classical spins was 
introduced 'by hand' such that the clusters are clearly formed. 
In this way, sites can easily be classified in three groups, as 
discussed below. Despite using a particularly chosen spin 
configuration, the resulting DOS is similar to those obtained 
previously (Fig|3j using a truly random distribution of spins. 



Lattice sites are classified as follows: the first group is defined 
to contain lattice sites that have a classical spin and belong 
to a spin cluster. The second class is composed of sites that 
have an isolated classical spin, i.e. nearby sites do not have 
other classical spins. Finally, lattice sites without classical 
spins in the same site or its vicinity belong to the third class. 
Ni(u), Nn(u), and N m (uj) denote the local DOS at sites 
corresponding to each of the three classes, respectively. It can 
be observed that Ni(u) contributes to states inside the impu- 
rity band (note the sharp peak near -1), but also contributes to 
the main band since the J/t used is intermediate. Nu(u) con- 
tributes a sharp peak near zero, and also has weight in the main 
band. This can be explained as follows: the electron is weakly 
localized at an isolated site yielding a state at ujh « — J or, 
since fi « — J for p « 0.4, ujh — fi « (see Fig|5k). This 
interesting result shows that many of the states near zero in 
Figs |3 and |4| may be localized, and they do not contribute to 
the conductivity. Finally, in Nm{uj) the empty sites contribute 
weight only in nearly unperturbed states, i.e., states that be- 
long to the main band. In conclusion, the local carrier density 
is very inhomogeneous and this fact has an important effect 
on the form of the site-integrated DOS. Our results show that 
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FIG. 3: N^(uj) and Ni(lj) vs. uj-ji for a 10x10 periodic system 
with 26 spins (x ~ 0.25) and 10 electrons (p ~ 0.4) at T/t = 0.01. 
Results are shown for (a) J/t = 2.0, (b) J/t = 3.0, (c) J/t = 4.0, 
and (d) J/t — 6.0. uj-\i is in units of the hopping, t. The results 
are averages over eight configurations of disorder (but there is no 
qualitative difference in the results from different configurations). 



(a) 



40 



20 



3 



-20 



-40 





" weak 


' 1 ' -J=2.0 


Up - 





\ Down 

V 

1/ 1 1 



3 -2 -1 1 2 3 4 5 
co-ji 



(c) 



40 



20 



3 



-20 



-40 



1 1 1 1 1 
. strong 


1 1 1 |-J=6.0 




— • — 
1 1 1 1 1 


Down^A - 
■ i i i i i i i i 



3 -2 -1 1 2 3 4 5 




-3 -2 -1 1 2 3 4 5 
co-ji 




1 2 3 4 5 6 
J/t 



FIG. 4: Nf(u) and Ni(u) vs. u-fi for a 6 3 periodic system with 
54 spins (x - 0.25) and 16 electrons (p - 0.4) at T/t = 0.01. 
Results are shown for (a) J/t = 2.0, (b) J/t = 4.0, (c) J/t = 6.0. 
uj-ji is in units of the hopping, t. The results are averages over eight 
configurations of disorder (but there is no qualitative difference in 
the results from different configurations), (d) Xc and T* vs. J/t 
for a 6 3 lattice with x — 0.25 and p — 0.4. These temperatures 
were determined from the spi n- sp in correlations at short and large 
distances, as explained in Ref.tl5l 
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FIG. 5: Local DOS for three different classes of sites (see text) on 
an 8x8 lattice with J/t = 2.5, T/t = 0.01, x ~ 0.2 and p = 0.4. 
This coupling regime corresponds to the intermediate optimal region, 
where the impurity band is not fully formed. "FM CLUSTER" is 
Ni(cu) in the text notation, while "ISOLATED SPIN" and "EMPTY 
SITE" are Nu(lj) and iYm(cj), respectively. 



FIG. 6: Coupling dependence of the conductivity at low temperature 
in two dimensions, (a) ct(cj) vs. to for a 10 x 10 periodic system with 
26 spins (x ~ 0.25), 10 electrons (p - 0.4), T/t = 0.01, and for 
different J/t's as shown, (b) Drude weight, D, vs. J/t for the same 
lattice and parameters as in (a). 



scanning tunneling microscopy (STM) experiments would be 
able to reveal the clustered structure proposed here, if it indeed 
exists, when applied to DMS materials. 



IV. OPTICAL CONDUCTIVITY 

In the previous section, it was argued that at the optimal 
coupling J/t - where Tc is maximized - the impurity band is 
not completely separated from the main band. This conclusion 
is supported by results obtained from the optical conductivity 
as well, which is shown in Fig. |6ja) for a two-dimensional 
lattice. This optical conductivity has two main features: (i) 
a zero-frequency or Drude peak and (ii) a finite-frequency 
broad peak which is believed to correspond to transitions from 
the impurity band to the main band, as argued below. From 
Fig. HJa), it is observed that inter-band transitions are not 
much relevant at weak coupling (J/t=1.0), but they appear 
with more weight at intermediate couplings (J/t=2.0 — 3.0), 
where Tq is optimal. It is worth remarking that at J/t strong 
enough, e.g. J/t=6.0, when the impurity band is well-formed, 
the finite-frequency peak is weaker in strength than for the op- 
timal J/t due to localization. In fact, if J/t were so large that 
the "impurity band" is completely separated from the unper- 
turbed band, then carrier localization would be so strong that 
inter-band transitions would not be possible. The Drude peak 
is plotted separately in Fig.|6jb) as a function of J/t at low 
temperatures. As expected, when J/t increases and localiza- 
tion sets in, the conductivity of the cluster decreases. Simi- 
lar results for the optical conductivity and Drude weight are 
found in three dimensions (see Fig. [J. 

The rise in absorption that appears for J/t > 2.5 at in- 
termediate frequencies is due to inter-band transitions, i.e. 
transitions from the impurity band to the unperturbed band. 
The frequency of this peak, uo inter, depends upon J/t as well 




FIG. 7: Coupling dependence of the conductivity at low temperature 
in three dimensions, (a) ct(cj) vs. to for a 6 3 periodic system with 
54 spins (x ~ 0.25), 16 electrons (p ~ 0.4), T/t = 0.01, and for 
different J/t's as shown, (b) D vs. J/t for the same lattice and 
parameters as in (a). 



as carrier concentration p, but for J/t > 2.5 (which is an 
intermediate value) and p = 0.4 (a realistic value), it was 
observed that 0.5t < uointer < l.Ot. The model param- 
eter t was previously estimated 15 to be t~0.3eV, yielding 
Q.lbeV < Winter < 0.3eV, in agreement with experiments 
(see below). 

The temperature dependence of cf(uj) is shown in Fig. [8j 
and it is as follows: for T < Tc, the Drude weight, D, de- 
creases in intensity as temperature increases and the same 
behavior is observed for the finite-frequency peak. When 
T > Tc, D increases slightly at first and then remains con- 
stant as the temperature is increased further. This slight raise 
in D is subtle and could be linked to the decrease in resistiv- 
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ity observed for T > Tq for intermediate J/t as explained in 
the next section, since the system is clustered just above Tq. 
Since both the optical conductivity and Drude weight decrease 
with temperature, to satisfy the sum rule the kinetic energy 
must also decrease as temperature increases. This is indeed 
observed (not shown) and is to be expected since localization 
of the wave function is stronger above Tq . 
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FIG. 8: Temperature dependence of the optical conductivity of model 
Eq. 0. Shown is <j(cu) vs. uj for a 10 x 10 periodic system with 26 
spins (x ~ 0.25), J/t = 2.5, p = 0.3, and for different tempera- 
tures, as indicated. Inset: Drude weight, D, vs. temperature T. 

A broad peak at around 0.2eV is experimentally observed in 
the optical conductivity of Gai-^Mn^As, as shown in Fig.|9] 
which is reproduced from Ref . y[ This feature has been ex- 
plained before in two different ways: (i) as produced by tran- 
sitions from the impurity band to the GaAs valence bandit 
or (ii) as caused by inter- valence band transitions. 6 Due to the 
frequency range and temperature behavior observed, our study 
supports the first possibility, i.e., that the finite-frequency peak 
observed for Gai-^Mn^As at around 0.2 eV is due to tran- 
sitions from the impurity band to the main band. Recent 
observations 34 also appear to support this view. 

The Drude weight increases in intensity with increasing 
carrier density, as can be seen from Fig. [TO] In fact, the ra- 
tio p of carriers to Mn concentration, which is a measure of 
the compensation of the samples, could explain, within the 
framework of this model, the different behavior observed at 
low frequencies for GaMnAs and InMnAs. In the former case, 
no tail of the Drude peak is found. However, for InMnAs, a 
clear Drude tailing is observedpi2& with increasing intensity 
as the temperature increases. At large enough doping, our cal- 
culations based on model Eq. predict a Drude-like peak 
for cf(uj) which could correspond to the regime valid for In- 
MnAs, while for low doping this study predicts a very small 
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FIG. 9: Absorption coefficient a(co) spectra for a metallic sam- 
ple prepared by low temperature annealing. The temperatures from 
down to top are 300, 250, 200, 150, 120, 100, 80, 60, 40, 20, 10 and 
4.2 K. (from Katsumoto et aim) 



Drude peak, which is the case for GaMnAs. Unfortunately, 
the precise carrier concentrations for both materials are not 
known experimentally with precision. 
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FIG. 10: p-dependence of the optical conductivity. Shown is 
including the Drude weight, vs. cu for a 10x10 lattice, x = 0.25, 
J/t — 2.5, T/t — 0.01, and different values of p. In the direction 
of the arrow, p takes the values 0.08, 0.12, 0.20, 0.30, 0.38, and 0.5. 
Inset: Drude weight, D, vs. fraction of carriers, p. 



V. CONDUCTANCE AND METAL-INSULATOR 
TRANSITION 

The conductance, G, is calculated using the Kubo formula 
adapted to geometries usually employed in the context of 
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FIG. 1 1 : Geometry used for the calculation of the conductance. The 
interacting region (cluster) is connected by ideal contacts to semi- 
infinite ideal leads. 



mesoscopic systems.- 55 The actual expression is 



G 



e 2 
2— Tr 
h 



(ihv x )Img(E)(ihv x )Img(E) 



(14) 



where v x is the velocity operator in the x direction and 
ImQ(E) is obtained from the advanced and retarded Green 
functions using 2ilmg(E) = Q R {E) - Q A (E), where E is 
the Fermi energy. The cluster is considered to be connected by 
ideal contacts to two semi-infinite ideal leads, as represented 
in Fig.[H] Current is induced by an infinitesimal voltage drop. 
This formalism avoids some of the problems associated with 
finite systems, such as the fact that the conductivity is given 
by a sum of Dirac S functions, since the matrix eigenvalues are 
discrete. A zero-frequency delta peak with finite weight cor- 
responds to an ideal metal - having zero resistance - unless 
an arbitrary width is given to that delta peak. In addition, in 
numerical studies sometimes it occurs that the weights of the 
zero-frequency delta peak are negative due to size effects^ 
For these reasons, calculations of d.c. resistivity using finite 
close systems are rare in the literature. All these problems are 
avoided with the formulation described here. 

The entire equilibrated cluster as obtained from the MC 
simulation is introduced in the geometry of Fig.[H] The ideal 
leads enter the formalism through self-energies at the left/right 
boundaries, as described in Ref.1351. In some cases a variant 
of the method explained in Ref . [35| was used in this paper to 
calculate the conductance (this modification was suggested to 
us by J. Verges). Instead of connecting the cluster to an ideal 
lead with equal hoppings, a replica of the cluster was used at 
the sides. This method, although slightly slower, takes into 
account all of the Monte Carlo data for the cluster, including 
the periodic boundary conditions. In addition, averages over 
the random Mn spin distributions are carried out. The phys- 
ical units of the conductance G in the numerical simulations 
are e 2 /h as can be inferred from Eq. (TRl . In Figure[l2|the in- 
verse of the conductance, which is a measure of the resistivity, 
is plotted for a three-dimensional lattice at weak, intermediate, 
and strong coupling, at fixed £=0.25 andp=0.3. For the weak 
coupling regime (J/t=1.0) the system is weakly metallic at 
all temperatures. In the other limit of strong couplings, \/G 
decreases with increasing temperature, indicating a clear insu- 
lating phase, as a result of the system being in a clustered state 
at the temperatures exploredpSi 5 . with carriers localized near 
the Mn spins. At the important intermediate couplings em- 
phasized in our effort, the system behaves like a dirty metal 
for T < T c , while for T c < T < T*,l/G slightly decreases 
with increasing temperature, indicating that a soft metal to in- 
sulator transition takes place near Tq. For T > T*, where the 
system is paramagnetic, 1/G is almost constant. Note that for 



strong enough J/t, Tq — > and therefore no metallic phase 
is present. 
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FIG. 12: Dependence of the theoretically calculated resistivity, p, 
with temperature in three dimensions. Shown is p = L/G vs. T 
on 4 3 lattices, 16 spins (x = 0.25), and 5 carriers (p — 0.3) for the 
J/t's indicated. An average over 20 disorder configurations has been 
performed in each case. Units are shown in two scales, ah/ (2e 2 ) on 
the left and mQcm on the right, with L — 4 and assuming a — 5.6 A. 
The estimated critical temperatures are also shown (the arrows indi- 
cate the current accuracy of the estimations). 

Similar qualitative behavior is found for the two- 
dimensional case (Figs. [T3l and [TH . Furthermore, in Fig.[T4l 
the spin- spin correlations have been plotted to show the loca- 
tion of Tc and T* and their relation to the resistivity. 

It is interesting to compare our results with experiments. 
For Gai-^Mn^As, data similar to those found in our investi- 
gations have been reported^ A typical result for resistivity 
versus T can be found in Ref. 4. The qualitative behavior of 
the resistivity in these samples agrees well with the theoretical 
results presented in Fig.[l2|if intermediate couplings are con- 
sidered. On the other hand, stronger or weaker couplings are 
not useful to reproduce the data, since the result is either in- 
sulating or metallic at all temperatures, respectively (defining 
metallic and insulator regimes by the slope of the resistivity 
vs. temperature curves). 

Furthermore, even the experimental numerical values of the 
resistivity are in agreement with the results presented here 
(even though model Eq. |l} does not include a realistic treat- 
ment of the GaAs host bands). This can be shown as fol- 
lows. The conductivity of a three-dimensional sample is re- 
lated to the conductance by a = G/L, where L is the side 
length of the lattice. Hence the resistivity is p = L/G. The 
units of G are, as explained before, h/(2e 2 ), and in our case 
L = 4a, where a is the lattice spacing. Assuming a = 5.6 
A, then the values shown in Fig. [12] are obtained. Note that 
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FIG. 13: Inverse of the conductance, 1/G, vs. T calculated on a 
lOx 10 lattice with 26 spins (x ~ 0.25), 8 electrons (p « 0.3), and 
two values of J/t as indicated. Shown is an average over three dis- 
order configurations. 1/G has units of h/(2e 2 ) in two dimensions. 



in Fig.[2|the minimum resistivity for J/t=2.5 is 3.3 m£7cm, 
whereas the minimum possible value for that carrier doping 
and cluster size used is 1.5 m^cm, which corresponds to the 
case J/t=0. In spite of the label 'metallic' for these results, 
the absolute values of the resistivity are high. Similarly, in the 
metallic phase of the sample shown in the experimental results 
of Ref. 4, the minimum resistivity is only ~ 3-6mOcm. Both 
in theory and experiments, the metallic phase appears to be 
'dirty', which is likely due to the reduced number of carriers, 
and localizing effect of the disorder. 

For J/t intermediate or strong, localization of the wave 
function is observed at intermediate temperatures^ This im- 
plies that carriers tunnel between impurity sites without vis- 
iting the main band, leading to insulating behavior. Below 
Tq conduction is favored by the ferromagnetic order and, 
as a consequence, as temperature decreases conductance in- 
creases. To further understand why the conductance (resis- 
tance) has a minimum (maximum) around Tq it is helpful to 
consider the three states: FM, clustered, and paramagnetic 
(PM) as depicted for a special spin configuration in Fig. [21 
on an 8 x 8 lattice. This configuration is not truly random but 
it was chosen 'by hand' for simplicity so that the spin clusters 
can easily be recognized. In the case of a random configura- 
tion of spins a similar reasoning can be drawn. For the FM 
state of Fig. [21 conduction is possible and indeed the mea- 
sured conductance is G w 2.8(2e 2 /h). For the two possible 
clustered states represented in Figs. [l3l b) and (c), the con- 
duction is reduced due to the different alignment of spins in 
different clusters. This is verified by calculating the conduc- 
tance which is G « 0.0 in those cases. For a paramagnetic 
spin arrangement with spin pointing in random directions (not 
shown), the conductance is small but finite, G w 0.5(2e 2 /h). 
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FIG. 14: Magnetization, \M\, spin-spin correlations, C(d), at max- 
imum and minimum distance, and inverse of the conductance, 1/G, 
vs. T on a 12x12 lattice with 22 spins (x ~ 0.15), 6 electrons 
(p ~ 0.3), and J/£=1.0. The approximate vanishing of spin- 
correlations at different temperatures, depending on the distance d 
studied, allowed us to obtain an approximate determination of Tc 
and T*, similarly as carried out in previous studies.— 

It follows that the clustered state has the minimum conduc- 
tance, and this explains the observed behavior of the resistivity 
with a maximum near Tq, i.e., when the system is clustered. 

The inverse of the Drude weight, D~ x , also shows a peak 
around Tq as seen in Fig. [OH where the same Mn-spin config- 
uration as in Fig. [14] was used. This provides further evidence 
for the metal-insulator transition near Tq described here. 

It is also interesting to remark that the optimal value of J/t 
changes with x (e.g., for small x, J/t optimal is also small). 
Since J/£ ptimai approximately separates the metal from the 
insulator, working at a fixed J/t (as in real materials) and 
varying x, then an insulator at high temperatures is found at 
small x, turning into a metal at larger x. This is in excel- 
lent agreement with the experimental results of Ref. 4 using a 
careful annealing procedure (while results of previous investi- 
gations with as grown samples had found an insulator-metal- 
insulator transition with increasing xm). The issue discussed 
in this paragraph is visually illustrated in Fig.[nl for the ben- 
efit of the reader. 



VI. MAGNETORESISTANCE 

The magnetoresistance percentage ratio was calculated us- 
ing the definition MR = 100 x (p(0) - p(B))/p(B). Fig- 
ure[n3 shows the magnetization, \M\, and the resistivity as a 
function of the applied field B on a 12x12 lattice, x=0.15, 
p=0A, J/t=1.0, and low temperature. Note that in two di- 
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FIG. 15: Qualitative explanation of the transport properties of clus- 
tered states. Shown are three clusters created 'by hand' on an 8x8 
lattice and with spin configurations also selected 'by hand' to illus- 
trate our ideas. For a (a) FM state at J/£=2.5 and p=0A, conduction 
is possible due to the alignment of magnetic moments and the con- 
ductance was found to be G=2.8(2e 2 /h). For the "clustered" state 
regime, two typical configurations (b) and (c) are shown where the 
conduction channels are broken and as a consequence G « 0. For 
the same clustered state but with randomly selected spin orientations 
(not shown) the conductance is small but finite G « 0.5(2e 2 //i). 
The radius of the solid circles represent the local charge density. 



mensions at £=0.25 andp=0.3, the intermediate coupling cor- 
responds to J/t=2.5, but at x=0.15, now J/t=1.0 corresponds 
to optimal coupling since as x decreases the optimal J/t also 
decreases. For these parameters, Tq was estimated to be 0.02£ 
and T* ^0.06£. The value in Teslas of the magnetic field 
was calculated assuming g=2.0 and 5=5/2 for the localized 
spins. The units of the resistivity are h/(2e 2 ) for the two- 
dimensional lattice. 

In our simulations it was observed that \M\ increases with 
increasing magnetic field. At zero magnetic field \M\ is 60% 
of its maximum value, while at 12 T it has reached ~80%. 
This is in agreement with the experimental results in Fig. 3 of 
Ref.0 where the magnetization is 50% of its maximum value 
at zero field, but it reaches 70% of that maximum value at 4 T. 
In our studies it was also observed that 1/G decreases with in- 
creasing magnetic field, showing at all fields a negative mag- 
netoresistance. Near Tq the resistivity decreases by 20%-30% 
increasing the field from zero to 12 T, while the decrease in re- 
sistivity at higher temperatures is much smaller. The present 
computational results agree very well with the experimentally 
measured dependence of the resistivity and magnetization on 
magnetic fields. For example, in Fig. 2(b) of Ref. 2, the de- 
crease in resistivity is 25% near Tq « 60 K when increasing 
the field from zero to 15 T. A much smaller decrease in resis- 
tivity is observed at higher temperatures (T > 120 K). It is 
also indicated in that experimental reference that the magne- 
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FIG. 16: Inverse of the Drude weight vs. T for a 12 x 12 lattice and 
parameters as in Fig. [l4| Note that both the inverse conductance in 
Fig-El an d Drude -1 provide evidence for a metal-insulator transi- 
tion at approximately the same temperature Tc . 
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FIG. 17: Dependence of J/£ op timai with x, approximately separat- 
ing the metal from the insulator (standard notation). At fixed J/t 
(horizontal solid line), as in experiments, the system should trans- 
form from an insulator to a metal with increasing x, in agreement 
with experiments 4 . 



toresistance is between 30% and 40% at 12 Tesla. 

It is interesting to remark that substantial magnetoresis- 
tance effects have been reported in thin films consisting of 
nanoscale MnnGeg ferromagnetic clusters embedded in a di- 
lute ferromagnetic semiconductor matrix. 37 The characteris- 
tics of these materials are analogous to our clustered state, that 
also show a robust MR effect due to magnetization rotation of 
spontaneously formed clusters. 



VII. MORE ABOUT CLUSTER FORMATION 

It was discussed before in this paper and in previous pub- 
lications that a clustered state is formed above Tc for inter- 
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FIG. 18: Magnetization, \M\, and inverse conductance, 1/G, vs. 
magnetic field, B, on a 12x12 lattice with PBC, x=0.15, p=0.3, 
J=1.0, and T/t=0.01 (same parameters as in Fig. ITU. The units 
of B are Tesla, assuming g=2.0 and 5=5/2. The units of 1/G are 
h/2e 2 . The approximate values of the magnetoresistance are 18% at 
4 Tesla, and 30% at 12 Tesla. 
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FIG. 19: Percolation picture for the formation of the clustered- state 
regime, (a) Two-dimensional lattice representing randomly-located 
classical spins as black dots with x=0.05. (b) Black areas represent 
non-zero carrier wave function, assuming a step-function profile for 
the wave functions with radius equal to 2 sites, as explained in the 
text, (c) Same as (b) but now showing connected regions (which 
could in practice correspond to the FM clusters discussed in the text) 
indicated with different shades of gray. 



mediate and large J/t couplings. This state is a candidate 
to describe DMS materials since it explains both the resistiv- 
ity maximum around Tq, as well as the decrease in resistivity 
with increasing applied magnetic field. In addition, it provides 
an optimal Tq . This clustered state is formed only for inter- 
mediate or larger J/t and when the compensation is strong, 
p < 0.4. Since for large J/t the carriers are localized, the 
problem becomes one of percolation theory, which has already 
been treated using different approaches and approximations. 24 
Here, only a very simple way of visualizing this state will 
be presented. First, consider a two-dimensional lattice with 
5% Mn spins represented as black dots, as shown in Fig.[T9h. 
The carrier wave function, ip(r), is considered to be localized 
around a Mn spin and it is assumed to be a step function for 
simplicity, i.e. ip(r) is non-zero only if r < r\ where r\ is 
the localization radius introduced by hand. Sites where the 
wave function is not zero are the black areas of Fig. fTW b). In 
this case, sites can be classified into connected regions, and 
that feature is indicated in Fig. [T9l c) using different shades. 
Each region is correlated and will correspond to a FM clus- 
ter. In this case all spins belong to some large cluster. As the 
concentration x grows, the clusters will tend to merge. As x 
decreases, these clusters will become more and more isolated. 

Recent experimental work on (Ga,Cr)As have revealed un- 
usual magnetic properties which were associated with the 
random magnetism of the alloy. The authors of Ref. 
explained their results using a distributed magnetic polaron 
model, that resembles the clustered- state ideas discussed here 
and in Refs. ll5l23l24l 



VIII. CONCLUSIONS 

In this paper, dynamical and transport properties of a single- 
band model for diluted magnetic semiconductors have been 
presented. The calculations were carried out on a lattice and 
using Monte Carlo techniques. The optical conductivity, den- 
sity of states, and resistance vs. temperature agree with ex- 
perimental data for Gai_ x Mn x As if the model is in a regime 
of intermediate J/t coupling. In this region, the carriers are 
neither totally localized at the Mn sites nor free, as in previ- 
ous theories. The state of relevance has some characteristics 
of a clustered state, in the sense that upon cooling from high 
temperatures first small regions are locally magnetized at a 
temperature scale T* - causing a mild insulating behavior - 
while at an even lower temperature Tq the alignment of the in- 
dividual cluster moments occurs - causing metallic behavior. 
This is in agreement with previous investigations ^i^SiSi 

One of the main conclusions of the paper is summarized in 
Fig|20| where the 'transport' phase diagram is sketched. The 
low-temperature ferromagnetic phase is metallic (although of- 
ten with poor metallicity). The clustered state between T* and 
Tq has insulating properties, and the same occurs in a good 
portion of the phase diagram above these two characteristic 
temperatures. Note that here the terms metallic and insulat- 
ing simply refer to the slope of the resistivity vs. temperature 
curves. The values of the resistivities in the two regimes are 
not very different, similarly as observed experimentally. 

This and related efforts lead to a possible picture of DMS 
materials where the inhomogeneities play an important role. 
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FIG. 20: Sketch of the phase diagram indicating the conducting and 
insulating regions, as obtained in the present investigations. The 
dashed area in the paramagnetic phase indicates the crossover re- 
gion from a mild metallic weak coupling regime to a mild insulating 
strong coupling regime. 



In this respect, these materials share characteristics with many 



other compounds such as manganites and cuprates, where cur- 
rent trends point toward the key importance of nanocluster 
formation to understand the colossal magnetoresistance and 
underdoped regions, respectively. "Clustered" states appear 
to form a new paradigm that is useful to understand the prop- 
erties of many interesting materials. 
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